Field Dependence of the Superconducting Basal Plane Anisotropy of TmNi2B2C 
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The superconductor TmNi2B2C possesses a significant four-fold basal plane anisotropy, leading 
to a square Vortex Lattice (VL) at intermediate fields. However, unlike other members of the 
borocarbide superconductors, the anisotropy in TmNi2B2C appears to decrease with increasing field, 
evident by a reentrance of the square VL phase. We have used Small Angle Neutron Scattering 
measurements of the VL to study the field dependence of the anisotropy. Our results provide a 
direct, quantitative measurement of the decreasing anisotropy. We attribute this reduction of the 
basal plane anisotropy to the strong Pauli paramagnetic effects observed in TmNi2B2C and the 
resulting expansion of vortex cores near i/c2. 

PACS numbers: 74.25.Uv,74.70.Dd,61.05.fg 
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I. INTRODUCTION 

The vortex lattice (VL) symmetry and orientation in 
clean type-II superconductors depends sensitively on the 
host material anisotropy, vortex density and tempera- 
ture, frequently leading to rich phase diagrams. As a 
result, VL studies can be used as a sensitive probe of the 
anisotropy of the superconducting state. 

In superconductors with a sufScient four-fold basal 
plane anisotropy, either due to the pairing symmetry or 
the Fermi velocity, the VL undergoes a successive se- 
ries of generic symmetry and orientational transitions as 
the vortex density is incrcascd>ii^ At low fields a dis- 
torted hexagonal VL is observed, oriented with the unit 
cell diagonal along a crystalline high symmetry direc- 
tion. As the field is increased the VL undergoes a first- 
order reorientation and symmetry transition to a rhombic 
phase with the unit cell diagonal rotated 45° with re- 
spect to the hexagonal VL phase. Finally, upon further 
increase of the field, the rhombic VL continuously trans- 
forms into a square symmetry. The transitions are driven 
by the growing importance of the four-fold anisotropy of 
the vortex-vortex interaction as the vortex density in- 
creases, explaining why further changes of the VL struc- 
ture are usually not observed once the square phase has 
been reached. There exists, however, two striking ex- 
ceptions to this behavior as seen in the superconductors 
TmNi2B2C and CeCoIns^"— In both of these materials, 
the VL undergoes the normal progression of symmetry 
transitions described above at low fields. However, the 
square phase VL is found to be reentrant, and the VL 
undergoes the same sequence of transitions but in the 
reverse order as the field is further increased. This in- 
dicates a reduction of the superconducting basal plane 



anisotropy in these materials at high fields and is the 
main objective of this report. 

In both TmNi2B2C and CcCoIns the superconduct- 
ing state is strongly affected by Pauli paramagnetic 
effects]^—"— Briefly, there is a significant spin-polarization 
of the unpaired quasi-particles in the vortex cores, result- 
ing in an increased amplitude of the modulation of the 
magnetic fieldjiSrJ^ With increasing field the vortex cores 
are also predicted to expand and become more isotropic, 
leading to the reverse sequence of VL transitions.— 

Here we present the results of small-angle neutron scat- 
tering (SANS) experiments to directly measure the evolu- 
tion of the basal plane anisotropy in the high-field square, 
rhombic and hexagonal VL phases. This is possible by 
measuring a large number of higher order VL reflections 
in a manner analogous to our previous study of non- 
magnetic LuNi2B2C (no Pauli paramagnetic effects and 
no reentrance of the square VL phase). -^^^ Our measure- 
ments allow a quantitative determination of the four-fold 
basal plane anisotropy and show a monotonic decrease 
with increasing field. 



II. EXPERIMENTAL 

The sample was a single crystal of TmNi2B2C of mass 
387 mg and dimensions 9.0 x 8.0 x 1.0 mm'^ with the c 
axis along the thin direction, grown by a high temper- 
ature flux method and using isotopically enriched ^^B 
to reduce neutron absorptionJ^ TmNi2B2C has a super- 
conducting critical temperature Tc = 11 K and a Neel 
temperature Tjy = 1.5 K below which the Tm moments 
order antiferromagnetically (AFM) J^"— The sample was 
mounted on an aluminium plate with the crystalline a 
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axis vertical. The c axis is horizontal and rotated by an 
angle f2 with respect to the applied magnetic field and 
the incoming neutrons as shown in the inset to Fig. [Ija) . 
The rotation (fJ) can favor a single VL domain orien- 
tation, thus reducing complications to the data analysis 
resulting from overlapping peaks from different domains 
while measuring higher orders of Bragg reflections. The 
SANS experiment was carried out at the 30 m NG7 in- 
strument at the NIST Center for Neutron Research, us- 
ing a neutron wavelength An = 0.55 nm and a spread 
AAn/An 11% FWHM. A horizontal field cryomagnet 
was used to reach the desired fields and temperatures. 
Measurements were done at a temperature of 1.6 K and 
in a field range 0.2 < ^qH < 0.6 T. Preliminary mea- 
surements were also carried out at the NG-2 SANS in- 
strument at Oak Ridge National Laboratory and at the 
Dll SANS at Institut Laue-Langcvin. 

Comparing the field dependence of the VL form factor 
(see Sect. Ill) for the first-order rcficctions to our previ- 
ous measurements showed a perfect agreement)^ confirm- 
ing that T > T/v as the AFM ordering of the Tm mo- 
ments significantly affects the VL form factor and sup- 
presses the Pauli paramagnetic effect;^ Vortex lattices 
were prepared by cooling through Tc in a constant field 
(FC) and, in some cases, followed by a damped small- 
amplitude field oscillations (FCO). Background mea- 
surements were measured at 14 K and were subtracted 
from the foreground measurements. The diffracted neu- 
trons were detected by a two-dimensional '^He position- 
sensitive proportional counter. 

III. RESULTS 

Examples of the square, and the high-field rhombic 
and hexagonal VL phases observed in TmNi2B2C are 
shown in Fig. [T] In this material the low-field hexago- 
nal and rhombic lattices occur below the convenient field 
range for SANS experiments. However, a hexagonal VL 
have been observed at 2 mT by Bitter dccoratioui^i^ For 
all three cases shown in Fig. [1] the measurements were 
extended to include as many higher-order reflections as 
possible within reasonable count times. The square VL 
diffraction pattern in panel (a) was obtained at 0.2 T 
and shows Bragg reflections with scattering vectors given 
by Qhk = {h^ + k'^Y'^qo, where go = 2^{B/^of'^ and 
(/)o = h/2e = 2068 T nm^ is the flux quantum. Fig. (Hd) 
shows the indexing of the peaks in one quadrant, with 
the remaining obtainable by symmetry. The diffraction 
pattern in Fig. [TJa) was obtained with O = 10°. Due 
to the very modest ac-anisotropy of the extrapolated or- 
bital upper critical field F = ^^(4'^oib/-^c2 orb ~ 1.2 in 
TmNi2B2C,— no distortion of the VL due to the rota- 
tion of the field away from the c axis is detectable. 

Increasing the field to 0.35 T causes the VL to change 
to a two-domain rhombic structure as shown in Fig. [IJb) 
with the corresponding indexing in panel (e) . The rhom- 
bic VL has an opening angle /3 = 73.0° and scatter- 



ing vectors q^k = {h"^ + k"^ + 2hkcos PY^"^ qo, where 
qo = 2TT{B/(f)osinf3)^''^. These measurements were done 
with O = 17°, chosen to favor one of the two VL domain 
orientations while at the same time keeping the distortion 
due to the ac-anisotropy small. For this the minority 
domain (red circles) is sufficiently suppressed to allow a 
reliable measurement of the intensity of the higher order 
majority domain (black circles) rcficctions. No measur- 
able difference in the VL opening angle (3 was observed 
between the two domains. 

Finally, as the field is increased to 0.5 T, a distorted 
hexagonal VL was observed as shown in Fig. [IJc) with 
the indexing in panel (f), and with an opening angle 
f3 ~ 56.1°. The magnitude of the scattering vector is 
given by the same expression as in the rhombic phase. 
Note that equivalent scattering vectors were chosen as 
the unit vectors, leading to qii being slightly shorter than 
qio as /3 < 60°. However, this is merely a naming con- 
vention and will not affect the analysis of the scattered 
intensity. The measurements were performed at a field 
rotation fl = 10°. For the hexagonal VL orientation the 
two domains are orientated equivalently with respect to 
the field rotation axis and are thus equally populated. A 
distorted hexagonal VL was also observed at 0.6 T (not 
shown) with an opening angle /S = 56.3°. However, at 
this field, no higher-order peaks were measurable due to 
the decreasing scattered intensity with increasing scat- 
tering vector and applied field. 

The diffraction patterns in Fig. [T] were obtained fol- 
lowing a preparation method chosen to produce the most 
ordered VL. Depending on the level and strength of vor- 
tex pinning in the host material relative to the vortex- 
vortex interaction, the optimal preparation may either 
be a field cooling (FC) procedure or a FC followed by a 
damped small-amplitude field oscillation (FCO) »^ While 
we found no difference between VLs prepared by the two 
different methods at 0.2 T, the FCO procedure provided 
a substantially better ordered VL at higher fields. This is 
seen in Fig. [21 which compares FC and FCO VL diffrac- 
tion patterns obtained at 0.5 T. It is likely that the VL 
disordering observed above 0.2 T in the FC case is due 
to the crossing of VL symmetry phase transitions while 
cooling from T > Tc to the measurement temperature of 
1.6 Ki^i^ Based on these findings, measurements were 
performed following a FC at 0.2 T and following a FCO 
for 0.35 T and above. In all cases where a FCO proce- 
dure was used, the initial amplitude of the damped field 
oscillation was 5% of the final field. 

Wc now turn to measurements of the VL form factors 
which are the main focus of this study. The form fac- 
tor F(qhk) is the Fourier transform, at scattering wave 
vector qhk, of the two-dimensional magnetic field mod- 
ulation due to the VL. It is related to the integrated 
reflectivity R, which is obtained by rotating and/or tilt- 
ing the cryomagnet and sample such that the VL scat- 
tering vectors cut through the Ewald sphere. Examples 
of rocking curves obtained in this fashion are shown in 
Fig. [31 In contrast to other members of the borocar- 



3 





-0.2 -0.1 0.0 0.1 0.2 -0.2 -0.1 0.0 0.1 0.2 -0.2 -0.1 0.0 0.1 0.2 

q^{nnT^) q^(nnT') q^{nm-^) 




FIG. 1. (Color online) SANS diffraction patterns showing the square, rhombic and hexagonal VL phases in TmNi2B2C at 1.6 K 
for applied fields of 0.2 T (a), 0.35 T (b) and 0.5 T (c) respectively. The images are obtained by summing measurements at 
multiple rotation and tilt angles to satisfy the Bragg condition for the different reflections. The scattered intensity is shown on 
a logarithmic false color scale to make strong and weak reflections simultaneously visible. The orientation of the crystalline axes 
and the magnetic field is shown in the inset to panel (a), where Q is the angle between the field and the c axis. The indexing of 
the VL Bragg peaks are shown for one quadrant in the schematics in panels (d) to (f) with the size of the circles indicating the 
intensity. For the rhombic and hexagonal VL phases two domain orientations with an opening angle /3 are observed, as shown 
by the black and red circles in (e) and (f). With increasing field the VL Bragg refiections move to longer scattering vectors q 
and decrease in intensity, making fewer of them visible. 



bide superconductors such as LuNi2B2C)iS the rocking 
curves in TniNi2B2C are found to be broad and with 
asymmetric line-shapes, necessitating multi-function fits 
to accurately obtain the integrated scattered intensity as 
shown for the (lO)-peak. For the higher-order reflections 
the line-shapes appear more regular, as seen for the (03)- 
pcak, and can be fitted by a single Lorcntzian, although 
this may also be a result of poorer signal-to-noise. Broad 
VL rocking curves, but with regular line shapes, were 
also found in other work, using TmNi2B2C single crystals 
from a different source<2^ It is possible that the present 
sample have more mosaicity which can explain the asym- 
metric lineshape. Nonetheless, the current rocking curves 
are still narrow enough to be easily measurable as shown 
in Fig. [3] and thus, the total scattered intensity can be 
precisely determined for each VL reflection. 

With the strong scattering from the VL in TmNi2B2C 
it is necessary to consider whether multiple scattering is 
affecting the measured intensities. This was discussed in 
detail by Densmore et al. in the case of LuNi2B2Cr^ 



where it was shown that multiple scattering did not pose 
a problem. In the case of TmNi2B2C, the integrated 
intensity is even stronger, but since the rocking curve 
is also significantly broader, the fraction of the incident 
neutrons scattered by the VL is < 0.4%, which is almost 
identical to LuNi2B2C. We thus conclude that the error 
in the measured intensities due to multiple scattering is 
insignificant. 

The integrated intensity is divided by the incident neu- 
tron flux to yield the integrated reflectivity 



Rhk = 



(1) 



where 7 — 1.913 is the magnetic moment of the neutron 
in nuclear magnetons and t is the average sample thick- 
ness (differences in t due to the change in Vt are < 3% and 
thus insignificant). The intensity for each reflection is 
corrected for the angle at which it cuts the Ewald sphere 
during the measurement of the rocking curve (Lorcntz 
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FIG. 2. (Color online) Comparison of VLs at 0.5 T {D, = 10°) 
and 1.6 K. The diffraction patterns were obtained directly fol- 
lowing a field-cooling procedure (FC) and after the applica- 
tion of a damped field oscillation with an initial amplitude of 
25 mT (FCO). No higher order VL reflections were observed 
due to shorter count times as compared to Fig. [TJc). 
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FIG. 3. Rocking curves for the square VL (10) and (03) Bragg 
reflections at 0.2 T and 1.6 K corresponding to Fig. [ija) and 
(d) . Note the different axes for the two reflections. Each angle 
was counted for 9 min. For the (10) reflection the error bars 
are smaller than the symbol size. The (10) reflection is fitted 
with a double Lorentzian function due to the irregular shape 
with a shoulder left of the main peak. The (03) reflections is 
fitted by a single Lorentzian. 



factor). Fig. 2] summarizes the measured VL form fac- 
tors for all reflections and fields of 0.2, 0.35 and 0.5 T. 
For the two higher fields the intensities of equivalent re- 
flections for the two domains have been added to obtain 
scattering from all the vortices in the sample in order to 
compare it directly to the square VL at 0.2 T. 




FIG. 4. (Color online) VL form factor divided by the applied 
field versus scattering vector q for all measured reflections at 
0.2, 0.35 and 0.5 T. For all fields and reflections the error bars 
are smaller than the symbol size. The curves are fits to the 
London model as described in the text. Full and dashed lines 
correspond to VL Bragg peaks along the crystalline [100] and 
[110] directions respectively. 



IV. DISCUSSION 



We will now discuss how the measured VL form factors 
can be used to study the evolution of the superconducting 
basal plane anisotropy in TmNi2B2C. Quahtatively, a re- 
duction of the anisotropy with increasing fleld is directly 
evident from Fig.|4l For the square VL at 0.2 T the form 
factors do not fall on a single curve, as expected for an 
isotropic superconductor, but rather lie on or between 
two curves going through \F{qhQ)\/H and \F{qhh)\/H, 
respectively. As these two limiting curves are 45° apart, 
their separation is a measure of the four-fold basal plane 
anisotropy. With increasing field we see that the separa- 
tions between the form factors along different directions 
decrease, indicating that the superconducting state be- 
comes more isotropic. However, the situation is compli- 
cated by the change in the VL symmetry, which changes 
the position of the reflections with respect to the crys- 
talline axes. 

The curves in Fig. 2] arc fits to the London model, 
extended by a Gaussian cutoff to take into account the 
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[110] 


[100] 




A (nm) ^ (nm) 


A (nm) ^ (nm) 


0.2 T 


64.1 6.28 


55.3 5.49 


0.35 T 




59.3 6.34 


0.5 T 


74.5 6.73 





TABLE I. Coefficients of London model fits (c = 0.5) shown 
in Fig. |4] in tiie case wliere the VL Bragg peaks are along the 
crystalline [100] or [110] directions. 



finite vortex core size^^i^ 

where A and ^ are, respectively, the penetration depth 
and the coherence length while c is a constant typically 
taken to be between 1/4 and Using c = 1/2, the 
fitted values of A and ^ for the curves in Fig. SI which 
correspond to crystalline high- symmetry directions, are 
given in Table HI It is important to stress that the primary 
objective of the fitting is to obtain an analytical expres- 
sion for the VL form factor for each field and direction, 
and the coefficients are not to be taken as accurate de- 
termination of X and ^. Nonetheless, Table |I] shows that 
increasing the field from 0.2 to 0.5 T both the fitted val- 
ues for the penetration depth and the coherence length 
along the crystalline [110] direction increases, consistent 
with a reduction of the field modulation and an expan- 
sion of the vortex core. The same is seen for A and ^ 
along [100] as the field is increased from 0.2 to 0.35 T. 

At 0.2 T, a measure of the superconducting basal plane 
anisotropy is obtained from the ratio (Aiio/Aioo)^ = 1.34. 
We note that a very similar value is found for the co- 
herence length ratio (^iio/Cioo)^ = 1-31. At the higher 
fields however, the changing VL symmetry causes the re- 
flections to move. Notably, at the two higher fields there 
are only VL Bragg peaks along one of the two crystalline 
high-symmetry directions ([100] for the rhombic VL at 
0.35 T; [110] for the hexagonal VL at 0.5 T). As a result 
it is not possible to directly extract the superconduct- 
ing basal plane anisotropy at the higher fields. In the 
following we present a more careful analysis of the field 
dependence of the anisotropy. 

As shown by Densmore et al.^ a conceptually simple 
and model-independent method to obtain a measure of 
the basal plane anisotropy is by a real space magnetic 
field reconstruction using 

B(r) = ^F(q„fe)e^'^-■^ (3) 

hk 

Since the SANS measurements only measure the abso- 
lute magnitude of the form factors, this requires an as- 
sumption about the relative sign of the Fourier com- 
ponents F{qhk)- In the case of LuNi2B2C, a compari- 
son to muon spin-rotation measurements showed that for 
fields below ^ ^fc2/3 the form factors all have the same 
sigUfi^ in agreement with the prediction of the London 
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FIG. 5. (Color online) Real space magnetic field reconstruc- 
tion from the measured VL form factors at /^o^^ = 0.2 T and 
T = 1.6 K. (a) Contour plot of the magnetic field, showing 
four VL unit cells with a vortex spacing of 98 nm, correspond- 
ing to a magnetic induction B — 0.216 T obtained from the 
magnitude of the scattering vector. Note that the image is 
rotated 45° with respect to Fig. [T] such that {110} directions 
are horizontal/vertical. The lowest contour corresponds to 
B = 198 mT and the contour spacing is 15 mT. (b) Current 
density as a function of distance from the vortex center along 
the VL nearest-neighbor direction ([110]) and the unit cell 
diagonal ([100]). The inset shows the value of £,j (distance 
of maximum current) in the basal plane. To emphasize the 
fourfold anisotropy, a circle with radius ^j^" is shown by the 
dashed line. 



model as well as numerical results based on the Eilen- 
berger equationsi^ In the case of TmNi2B2C the mea- 
sured Hc2 = 0.75 T at 1.6 K is severely Pauli limited<^ 
Instead, we use the extrapolated orbital upper critical 
field -ffc2.oib = 4.3 Tf2^ yielding an estimated upper limit 
of 1.4 T for the all-equal sign scheme. This is thus ex- 
pected to be valid for all the measurements in this re- 
port. Fig. [5]Ja) shows the real space field reconstruction 
obtained from the measured form factors in an applied 
field of 0.2 T. The accuracy of the reconstruction de- 
pends on the number of reflections included in the sum 
in Eq. ([3]). In the present case, the magnitudes of the 
form factors for the (30), (31) and (32)-rcflections, which 
were the highest order peaks which could be measured. 
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are less than 1 mT and do not change the reeonstruction 
to any significant degree. We note that the vortex spac- 
ing d = 27r/gio can be determined from the magnitude 
of the VL scattering vector, which aUows for a determi- 
nation of the magnetic induction B = <j>o/(P = 0.216 T 
in good agreement with earher reports.— That B > ^qH 
is due to the paramagnetism of TmNi2B2C for T > T/v. 
The magnitude of the field modulation w 130 mT or 0.65 
is much larger than the 10% observed in non-magnetic 
LuNi2B2C in an applied field of 0.5 T. Extrapolating 
the LuNi2B2C form factors to a field of 0.266 T yields 
a modest estimate for the increase of the field modula- 
tion: exp[27r2(8.22 nm)^ (0.5 T - 0.266 T)/^o] = 1.16-i^ 
This is still much smaller than the 60% observed here 
for TmNi2B2C. This difference is a manifestation of the 
strong Pauli paramagnetic effects, which leads to a sig- 
nificant polarization of the unpaired quasi-particle spins 
in the vortex cores of TmNi2B2C and creates a periodic 
magnetization that adds significantly to the field modu- 
lation from the circulating supercurrents. 

From the field reconstruction, one can calculate the 
current distribution hy ~ V x B, which con- 

tains contributions from both the supercurrents and the 
periodic magnetization and which can not be easily 
deconvoluted>i2iH Fig. [SJ^b) shows \ J[r)\ as a function 
of the distance from the vortex core along the VL near- 
est neighbor direction ([100]) and the unit cell diagonal 
([110]). The distance from the vortex center to the peak 
of the current density provides a measure of the core 
size Or^i^ which is seen to differ for the two direc- 
tions shown. The inset to Fig. EJb) shows the vortex 
core size in the plane perpendicular to the field, which 
displays a clear four-fold anisotropy. Specifically, we find 
= 12.4 nm and = 11.7 nm. The ratio between 
these two values is 1.06, slightly smaller than the 1.08 
found for the square VL in LuNi2B2C.— Since the VL at 
0.2 T is close to the onset of the square-to-rhombic tran- 
sition in TmNi2B2C, the measured anisotropy provides 
an estimate of the critical value necessary for stabilizing 
a square symmetry. Comparing the values for to the 
London model fits at 0.2 T listed in Table U one finds 
a substantial difference, unlike our earlier measurements 
on LuNi2B2C where the two were found to be in excellent 
agreement A theoretical analysis of the interdepen- 
dence of the superconducting and magnetic properties of 
TmNi2B2C by Jensen and Hedcgard gives an estimate of 
the orbital critical field i?c2,orb = 4.3 T at 1.6 yield- 
ing a zero-field coherence length of 8.75 nm. However, 
this ignores the contribution to the core size from the 
spin-polarized quasi-particles in the vortex core and we 
thus expect that obtained from the field reconstruction 
to be a more accurate measure of the actual vortex core 
size. It should also be noted that the anisotropy of O 
is smaller than the ratio from Table HI again illustrating 
that while ^ and are related they are not identical. 

While a field reconstruction at higher fields would be 
desirable, it is not possible to measure enough higher- 
order reflections necessary to obtain this with satisfactory 



accuracy. This is the case particularly in the hexagonal 
VL phase. Instead we return to the London model fits in 
Fig. [4] to obtain a quantitative, symmetry independent 
measure of the field dependence of the superconducting 
basal plane anisotropy in TmNi2B2C. The fits indicate 
that the VL form factor along any direction in the basal 
plane may be parameterized by 



H 



1 



(4) 



where (j) is the angle relative to the crystalline [110]- 
direction. This yields an expression for the 0-dependence 
of A and thus an angle dependent anisotropy ratio 



/AnoV 



e-c[g?(0)]2 _F(^)- 
This ratio is expected to exhibit a four-fold symmetry 

1 — cos 4fl!) 



(5) 
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(6) 



where the parameter a is the anisotropy amplitude. As 
noted earlier, the fitted values for A and ^ given in Ta- 
ble U show (Aiio/Aioo)^ ~ (^iio/6oo)^, indicating that 
A((/)) and have the same anisotropy amplitude. From 
Eq. ([U we obtain 
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In Fig. [6] we show the anisotropy ratio obtained by Eq. ([5]) 
using the measured VL form factors at 0.2 T and ^(0) 
calculated using Eq. (O and the values in Table HI The 
data are well fitted by Eq. ^ yielding an anisotropy 
amplitude a = 0.368 ± 0.014 in good agreement with the 
result based solely on the form factors corresponding the 
VL Bragg peaks on the high symmetry [110] and [100] 
directions (Table |T|. 

Having demonstrated that the approach above yields 
a consistent results for the anisotropy we now apply it 
to the 0.35 and 0.5 T measurements. Here we simul- 
taneously adjust the value of a used to calculate ^{4>) in 
Eq. dSI and the fitted value obtained by Eq. (|6]) to ensure 
a self consistent result. This approach allows a determi- 
nation of the anisotropy amplitude even in cases where 
there are only VL Bragg reflections along a single crys- 
talline high symmetry direction ([110] or [100]). Figure[6] 
shows the results of this analysis, providing a quantitative 
measure of the monotonically decreasing superconduct- 
ing basal plane anisotropy with increasing field. This is 
in stark contrast to the behavior found for non-magnetic 
LuNi2B2C, where a similar analysis on the data from 
from Ref. [la yields a = 0.456 ± 0.060 at 0.5 T increasing 
slightly to 0.492 ± 0.008 at 1.0 T. 
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FIG. 6. (Color online) (a) Angle dependence of the super- 
conducting basal plane anisotropy ratio calculated from the 
VL form factors as described in the text. For each field the 
curves show a fit to the anisotropy function given in Eq. ((6|. 
(b) The same results in polar coordinates. 

Our results provide a direct confirmation of theoreti- 
cal predictions that in superconductors with strong Pauli 
paramagnetic effects, paramagnetic dcpairing causes the 
vortex cores to expand and also become more isotropic 
as one approach H^i^^ One would therefore also ex- 
pect that the anisotropy of TmNi2B2C will increase be- 
low the antiferromagnetic ordering at T/v = 1.5 K where 



the Pauli paramagnetic effects are known to decreasei^ 



V. CONCLUSION 

In summary, we have investigated the field dependence 
of the superconducting four-fold basal plane anisotropy 
of TmNi2B2C. We have observed and quantified the de- 
creasing anisotropy with increasing applied field, which 
provides an explanation of the reentrant square VL 
phase. The decreasing anisotropy is attributed due 
to the strong Pauli paramagnetic effects observed in 
TmNi2B2C, leading to an expansion of the vortex cores 
near iJc2- We believe that a similar mechanism is respon- 
sible for the reentrance of the square VL phase observed 
in CeCoIns. 
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